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Abstract 

If the space Q of quadratic forms in R n is splitted in a direct sum Qi (B ■ ■ ■ ® Qk and if X 
and Y are independent random variables of R n , assume that there exist a real number a such that 
E{X\X+Y) = a(X+Y) and real distinct numbers 61, b k such that E(q(X)\X + Y) = b iq (X + Y) for 
any q in Qj. We prove that this happens only when k = 2, when R" can be structured in a Euclidean 
Jordan algebra and when X and Y have Wishart distributions corresponding to this structure. 



I Introduction 

Let S r be the set of (r, r) real symmetric matrices and let X and Y be independent random variables 
valued in S r such that they are Wishart distributed j p . a and 7 P ', CT , which means that 

E( e - tiex )=det(I r + 6a)- p (1.1) 

where 9 and a are in the set P r of the positive definite elements of S r and p is in 

.1 r — 1, ,r — 1 , 
A = {-,...,— }U(— ,00) (1.2) 

(In (ll.ip tr means trace). Note that for a = p/{p + p') 

E{X\X + Y) = a{X + Y). (1.3) 

Assume furthermore that p + p' > This implies that (X + Y)^ 1 exists. Then it is known that 

Z = (X + r)" 1 / 2 A(A + Y)- 1 ' 2 and X + Y are independent and that Z - uZu T for any orthogonal 
(r, r) matrix u. There are many consequences, nuances and characterizations of the Wishart distributions 
related to this result. One of these consequences is the following fact: for any s € S r consider the two 
quadratic forms on S r defined by 

<7i(x) = — tr 2 (xs) + tr (sxsx), q^{x) = tr 2 (xs) — tr (sxsx) (1-4) 

and the two numbers 

p p+l t p p- \ 



fo i = —, , r— r . & 2 = 



p + p'p + p' + l' p + + 
Then for { = 1,2 and for any s 



E(rf(X)|X + Y) = &ig?(X + y) (1.5) 
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This is the particular case d = 1 of Corollary 2.3 of Letac and Massam (1998). An important fact about 
this set (gfjC^seSr i s that it spans the whole space of quadratic forms Q on S r (since if q s (x) = tr 2 (xs) 
then {q s ; s £ S r } spans Q). More specifically denote by Qi the subspace of Q generated by {qf ; s € S r }- 
Then Q = Qi © Qi (see for instance Theorem 5.2 below for a proof). 

The aim of the paper is to prove a reciprocal statement of (|1.3p and (|1 .5|) : Let V be a linear real 
finite dimensional space (instead of S r ) and denote by Q the space of all quadratic forms on V. Fix 
a decomposition Q = Qi © Q 2 © ■ • ■ © Sfe with fc > 2 as a direct sum of linear subspaces. Consider 
two independent random variables X and Y with exponential moments satisfying (|1.3|) for some a and 
E(g(X)|X + Y) = biq(X + Y) for all q & Qi and for some distinct real numbers , &&. We show 

that under these circumstances, necessarily fc = 2 and X and Y are Wishart distributed in the following 
sense: there necessarily exists a structure of Euclidean Jordan algebra on V (like symmetric matrices, 
Hermitian matrices, or space with a Lorentz cone) such that X and Y are Wishart on the symmetric 
cone associated to it. Section 5 contains more detailed information about the two spaces Q\ and Qi of 
quadratic forms on S r (or more generally, on a Euclidean Jordan algebra) 

II Some history of the subject 

Wishart distributions on S r . Wishart distributions have been introduced by J. Wishart (1928) as 
distributions of Z\Z^ + • • • + ZjyZj^ ~ 7jv/2,2£ where Z\, . . . , Zjss are iid in W such that Zi ~ N(0, E). 
Elegant calculations about them are in Bartlett (1933) and the classical reference is Muirhead (1982). For 
the space S r of (r, r) real symmetric matrices the extension of the definition of 7 PiCT from a half integer p 
to the whole set A defined by (jl.2p is made in the fundamental paper of Olkin and Rubin (1962). Proving 
that a distribution 7 P!(T on the semi positive definite matrices such that (|1.1[) holds only if p is in A was 
considered as a challenge by statisticians (see Eaton (1983)) although the appendix of Olkin and Rubin 
contains an unnoticed proof of it (and unfortunately erroneous: see Casalis and Letac (1994)). This 
conjecture was independently proved by Shanbhag (1988) and Peddada and Richards (1989) by quite 
different means, although a solution already appeared in Gyndikin (1975) and seems to have been well 
known by analysts, who also call the set A and its extensions the Wallach set (see Lassalle (1987) for 
proofs and references). 

Lukacs-Olkin-Rubin Theorem. Wishart distributions on S r are the most natural generalization 
of the gamma distributions on the positive line. Lukacs (1956) shows that if X and Y are positive, 
independent non Dirac random variables and if Z = Xj(X + Y), then Z and X + Y are independent 
if and only if there exists u,p,p' > such that X ~ j p>a - and Y ~ 7p',o-- This was extended to S r by 
Olkin and Rubin (1962) by a proper definition of Z such that Z is symmetric (for instance by choosing 
Z = (X + Y)-!/ 2 A(A + Y)- 1 ' 2 or by choosing Z = C- X X(C- 1 ) T where C is the triangular matrix with 
positive diagonal elements coming from the Cholesky decomposition CC T = X + Y). They show that if 
X and Y are independent non Dirac random semi positive definite matrices in S r such that X + Y is 
invertible and such that Z ~ uZu T for any orthogonal (r, r) matrix u then Z and X + Y are independent 
if and only if there exists a positive definite matrix a and p and p' in A with p + p' > (r — l)/2 such that 
X - 7 p !(T and Y - 7 pV . If Z is defined as (X + Y)" 1 / 2 A(A + Y)^ 1 / 2 , Bobecka and Wesolowski (2002) 
have shown that the invariance hypothesis for Z by the orthogonal group can be dropped provided one 
assumes that X and Y have smooth densities. Removing this assumption of density is still a challenge. 

Wishart distributions on Hermitian matrices and on Euclidean Jordan algebras. Since 
normal distributions on Hermitian spaces have been considered (see e.g. Goodman (1963)), therefore 
Wishart distributions on Hermitian matrices occur naturally. Actually physicists considered them quite 
early (see Mehta (2004)). Carter (1975) in an unpublished PhD thesis extends Olkin and Rubin to this 
case. 

On the other hand, works on the classification of natural exponential families by their variance function 
have led to the observation that the exponential family {7 P , CT ; o € P r } of Wishart distributions on S r with 
fixed shape parameter p € A has a variance function which is the map from S r into itself x 1— > V(m)(x) = 
-mxm where m is in P r . In other terms, this means that if k is a cumulant function of 7 P . (J then for all 
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x in S r we have 



k"(0)(x) = -k'(6)xk'(6). 



Facts about multivariate distributions such that their corresponding variance functions are quadratic 
in the mean are collected in Letac (1989). In particular, Wishart distributions obtained from simple 
Euclidean Jordan algebras are described there. An indispensable reference for simple Euclidean Jordan 
algebras is Faraut and Koranyi (1994) always abreviated F.-K. below. Recall that simple Euclidean 
Jordan algebras are basically in one to one correspondence with the irreducible symmetric cones (self 
dual cones in Euclidean space such that the group of automorphisms of the cone acts transitively on it), 
in the way that S r is linked to P r . A quick definition of the Wishart distribution 7 P!<T on the Jordan 
algebra V with rank r, Peirce constant d, cone f2 of square elements, trace and determinant function tr 
and det can be done by its Laplace transform 



where a is in the interior £1 of il and where p is in the Gyndikin set of the Jordan algebra V defined by 



While the definition of determinant is the standard one for S r and for Hermitian matrices, it requires 
some care for the three other types of Jordan algebras: quaternionic Hermitian matrices, 27 dimensional 
Albert algebra and the algebra of the Lorentz cone. 

Particular cases of use of Wishart distributions on Jordan algebras in statistics occurred earlier (An- 
dersson (1975) for the Hermitian and quaternionic cases, and Jensen (1988) for the Lorentz cone, with 
its deep connexions to Clifford algebras). Jordan algebras are the natural framework for Wishart dis- 
tributions: Casalis and Letac (1996) is a clarification and an extension to Jordan algebras of Olkin and 
Rubin (1962) and of Carter (1975); Carter follows step by step the difficult Olkin and Rubin's approach 
and his work was unknown to Casalis and Letac (1996). 

Quadratic homogeneity and Wishart distributions. A remarkable fact about the classical 
Wishart distributions on S r is that the above variance function m V(m) is not only quadratic in 
m but homogeneous quadratic. This happens also to be true for Wishart distributions on any Euclidean 
Jordan algebra. This observation lead Casalis (1991) to prove the converse: any natural exponential 
family with a homogeneous quadratic variance function is a Wishart family, as conjectured in Letac 
(1989). Put in other words, if k is a cumulant function of some random variable X valued in M™ such 
that k"(9) = V(k'(8)) where V is a homogeneous quadratic function, then R" can be structured in a 
Jordan algebra such that X is Wishart for that structure. 

Quadratic REGRESSION property. A slight extension of Lukacs (1956) is to take two non Dirac 
independent rv X and Y on the positive line such that there exist positive a and b such that E,(X\X +Y) — 
a(X + Y) and E(X 2 \ X + Y) = b(X + Y) 2 and to prove that there exist positive p,p',cr such that X ~ 7 P)CT 
and Y ~ j p ,^. To see this, just multiply these two equalities by e e ^ x+Y \ take expectations and obtain 
two differential equations for the Laplace transforms of X and Y. This procedure is contained in Laha 
and Lukacs (1960). Bivariate regression version of Lukacs theorem based on conditions E(Xf\X + Y) = 
b(Xi + Yi) 2 , i = 1,2, where A = (Ai,A2) and Y = (Yi,!^) are independent was obtained in Wang (1981). 
This result was generalized in Letac and Wesolowski (2008) by considering regressions of quadratic forms 
E(q(X)\X + Y) = bq(X + Y) for all quadratic forms q orthogonal to an arbitrary fixed quadratic form 
go- That is in the setting of the present paper we required k — 1 and codimension of Qi to be equal 1. 

Letac and Massam (1998) use the quadratic regression approach to get a simpler proof of Olkin and 
Rubin theorem, as extended to Jordan algebras in Casalis and Letac (1996). It actually characterizes the 
Wishart distributions of independent X and Y in S r (and more generally of a Jordan algebra) through 
the following properties: if for i = 1,2, s G S r and qf are defined by (|1.4[) . then p.5[) holds (with 
suitable analogues of qi if the Jordan algebra is not S r ). Note that this regression perspective leads to a 




Ay = {^d,...,~(r-l)}U^(r-l),oo). 



(2.6) 
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characterization of 7 P , CT , 7p',<t without the hypothesis of invertibility of X + Y which was needed in the 
Olkin and Rubin characterization. 



Ill Main result 

Let V be a real linear space with dimension n > 1, let V* be its dual and consider the space T = 
L S (V,V*) of the symmetric linear maps from V to V* . If 9 £ V* and x £ V we write (6,x) for 9{x). 
Denote by Q the space of quadratic forms q on V, namely the set of real functions q on V such that 
(x, y) i — y \{q{x + y) — ?(x) — q(y)) is bilinear on V x V and q(Xx) — X 2 q(x) when A is a real number. The 
map from J- to Q defined by / t— > qf where x i— > Qf{x) = (f(x), x) is one to one. More specifically: 

\{q f {x +y)- q f (x) - q f (y)) = ~«/(x), y) + (f(y),x)) = (/(*), y) 

For q £ Q we therefore define the inverse map q i-> f q of / i-> qf by 

^(q(x + y) - q(x) - q{y)) = (f q {x),y). 

Let us also define here the concept of irreducibility for a probability measure /i on V. We say first 
that n is reducible if there exists a direct sum Vt ® V2 = V with diml^ > for i = 1,2, two probability 
measures [i\ and 112 on Vi and V2 such that /i = /ii (g) /i2- In other terms, if X ~ fi its projections Xi on 
V\ parallel to V2 and X2 on V2 parallel to V\ are independent. Suppose that furthermore X has a Laplace 
transform L — e K defined on some open set C V* = Vy ® V 2 * ■ In this case n(6) = Ki(^i) + ^2(^2) 
where Oi is the projection of 9 on V* and m and K2 are the cumulant functions of X\ and Xi- We say 
also that X and k are reducible in that case. Finally, /i, X and k are said to be irreducible if they are 
not reducible... 



Theorem 3.1 Let Qi Q2 © ■ • • © Sfe = Q be a direct sum decomposition of the space of quadratic 
forms on V with k > 2. Let X and Y be two independent irreducible random variables valued in V such 
that their Laplace transforms exist on an open set 9 c V*. We assume that 

1. there exists a real number a such that E(X|X + Y) = a(X + Y); 

2. there exist distinct numbers b\, . . . , bk such that for any i = 1, . . . , k and for any q £ Qi we have 

E(q(X)\X + Y) = b t q(X + Y). (3.7) 



Under these circumstances 0<a<l, k = 2 and there exists a simple Euclidean Jordan algebra structure 
on V such that X and Y are Wishart distributed on the positive cone of the algebra with the same scale 
parameter and respective shape parameters p and p' in Ay defined in (|2.6[) . Moreover Qi and Q2 are 
spanned by 

qt(x) = -tr 2 (xs)+tr(¥(x)(s)s), g|(x) = tr 2 (a;s) - tr (P(x)(s)s) (3.8) 

where tr , P and d are respectively the trace, the quadratic map and the Peirce constant of the Jordan 
algebra and s £ V. In this case 

P t P P+l t p P~ I 



p + p' ' P + p' P + p' + 1 ' P + P'p + p' — 5 



Proof. Denote by Lx and the Laplace transforms of X and Y. It is standard to prove that from 
condition 1) we have L]f a = L y : just multiply both sides of E(X\X + Y) = a(X + Y) by e< 9 > x + r > where 
€ O and take expectations of both sides to obtain the differential equation aL' x /Lx = (1 — a)L' Y / Ly ■ 
The fact that X and Y are irreducible implies that a = or a — 1 is impossible. The fact that log Lx and 
log Ly are convex implies that a < or a > 1 are impossible. From now on we denote e K — Lx — Ly^ 1 ^ ■ 
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In the sequel, we use the symbol Tr for the trace of an endomorphism. The symbol tr is reserved for 
the trace in a Jordan algebra. If q is a quadratic form on V we write 

q(^M(9) = Tr(f qK "(9)). 

Since k is a real twice differentiable function defined on an open subset of V* , the second derivative 
k"(9) is an element of L S (V* , V), the linear map f q is an element of L S (V, V*) and thus f q n"{ff) belongs 
to L(V* ,V*). It therefore makes sense to speak of the trace of this endomorphism of V* . Note that 
(f q (x),x) = Tr (f q (x £g> x)) and that <g> -gaK — n" . This explains the definition q{-§g){n) = Tr (f q K,"). 
Also q(K') can be written in terms of f q as <?(«/) = Tr {f q {n' ® «')) = k'). 

Calculations done in Letac and Wesolowski (2008) (2.9), show that for any i — 1, . . . k and for all 
q G Qi we have 

(i--)9(4)(«) = (^-i)g(«')- (3-io) 

(Again, to prove (|3.10p just multiply (|3.7[) by e ! < e ' X+Y ' > and take expectations). Observe that bi = a is 
impossible, since it implies that q(n') = for any q in Qi. Since Qi is not the zero space, there exists a 
non zero q with q(n') = 0. Now {x G V^; q(x) = 0} is a quadric of V and has an empty interior. On the 
other hand, since X is irreducible, this implies that X cannot be concentrated on some affine subspace of 
V. Therefore k is strictly convex and the set k'(@) is open and cannot be contained in a quadric. Thus 
a = bi is impossible, division by (1 — ^) is permitted and we rewrite (|3.10[) as 

?(^)(«) =**«(*') (3-H) 

where Vl = 

Now let us fix 9 G ^* and consider the element 9 ®9 oi F defined by (9 ® 0)O#) = (9, x)9. Denote by 
Ti the image of Qi by the isomorphism q i— > / g . Obviously we have 

J 7 ! © J" 2 e . . . © Jit = J 7 . 

Therefore there exist elements fa G J"j such that f\ + . . . + fu = 9 ® 9. Since /i, . . . , depend actually 
on we rather write /i(#, x) instead of fi(x) for x G V. Thus a; i— > /i(#, x) is a linear map from to V* . 
We rewrite the equality 9 ® 9 = f\ + . . . + fk as 

(M) 2 = (/x(M),x) + . . . + (/ fc (M), x) 

for any x in V. We now fix 9 — 9q in this equality and we recall that q(-§g)(n)(9) means Tr (f q n"(9)). 
Thus we apply the equality 



to k, the log of L. We get 



Tr ((9 ® fl ) = ^ Tr [/ 4 (# , •)«"(*)] • 



We now use the fact that x <-> (fi(9a, x),x) — q(x) is a quadratic form belonging to Qi to which we apply 
(pTTTj) . We therefore get 

Tr((0 o ®#o) K"(d))=X)piT>(/i(do,-)K / (fl)(»/5 / (fl)). (3.12) 
i=l 

Since this is true for any 9q in V* this is enough to claim that n" is a quadratic homogeneous function 
of k'. 
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We now apply the Casalis' theorem (1991) which says that if k is irreducible and if k" is a quadratic 
homogeneous function of k' , then there exists a simple Euclidean Jordan algebra structure on V related 
to X in a way that we explain now. Let O be the open cone of the squares of V, let tr and det be the 
trace and determinant functions on the Jordan algebra, let d and r be the Peirce and rank constants of 
V. Then there exists p <E Ay defined by (|2.6p and a € fi such X has the Wishart distribution j Py(7 on £1 
denned by its Laplace transform E(e" tr6x ) = det(J r + 0cr)- p for all 6 € O. 

To complete the proof, denote for a while by Q\ and Q 2 the spaces of quadratic forms spanned by 
(ql) s& v and (q 2 ) s ev as defined in (|3.8[) . Denote also 

~ P p+1 ~ p p- I 

01 = , 62 = ^-j- 

p + p' p + p' + 1 p + p' p + p' - I 

Recall that we want to prove that k = 2 and that {Qi, Q2} = {Qi, Q2j-- Let now g 6 Qi. Therefore 
E(q(X)\X + Y) = biq(X + Y). We now write q — qi +92 with e Qi which is possible since Qi© Q2 = Q- 
Recall that since X and Y have distributions 7 p cr and 7 P ' lC r we can write E{qi{X)\X + Y) = b L q(X + Y). 
Thus 

(61 - h) qi {X + Y) = {h - b 2 )q 2 (X + Y). 

Since X + Y is valued in the open set fi this implies (b\ — bi)q\ = {bi — b 2 )q2- Thus the two sides of this 
equality are zero: either bi — b\ and q 2 — or the reverse statement holds. Since we have assumed that 
bi, . . . , bk are distinct, this ends the proof. 



IV Comments 

1. Surprisingly enough, while starting from a linear space V without any additional algebraic structure, 
the regression conditions on X and Y of the theorem impose by themselves a Euclidean Jordan 
algebra structure on V. 

2. The three numbers a, 61 and b 2 together with the dimension of V determine uniquely the structure 
of Jordan algebra on V in the following sense: we can see from the equations (|3.9[) that b 2 < a 2 < 
bi < a. Moreover these equations give the Peirce constant d of V by 

d=2 a-b! a 2 -b 2 



b\ — a' 2 a — b 2 

Since the rank r satisfies dim V = r + ^r(r — 1) the type of the Jordan algebra is completely known. 

In the theorem, k = 1 would lead to X and Y concentrated on a line M.v of V. If X = X\v 
and Y = Y\v then X\ and Y\ would be one dimensional gamma distributed and X would not be 
irreduciblesince we have assumed diml^ > 1. Furthermore if in the theorem we do not assume that 
&i, . . . , bk are distinct, then either they are all equal to one b and this sends us back to the trivial 
case k = 1 or they are not and if k 1 > 2 is the number of distinct b^s, then the theorem gives k' = 2. 

Some comments about irreducibility are in order. If Ly is a power of Lx, then Y is irreducible if 
and only if X is. Therefore irreducibility can be assumed in the theorem for X only. If irreducibility 
is not assumed, we have an artificial generality. For instance suppose that {Xi,X2,X^, Y^Y^^a) 
are independent real rv with Xi ~ 7 Qi ,o- and Yi ~ 7ft ,<7- Then for i 7^ j we have 

E(X lXj \X + Y) = ^2 -n-A Xi + YXXj+Y,), 

{at + Pi){aj + fij) 

E{X 2 \X + Y) = ^'\T\ 1 W^ + Y ^- 
{a, + Pi) [pa +Pi + 1) 

This implies that k — 6 corresponding to the 6 independent quadratic forms on V — R 3 defined by 
qij{x) = XiXj for i < j. 
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V The spaces Qi and Q2'. the operator \I/ 



If V is a simple Euclidean Jordan algebra with rank r and Peirce constant d — 2d', denote by T — L S (V) 
the space of symmetric linear operators on V. The dimension of V is n = r + dr(r — l)/2. Given y £ V, 
important examples of elements of T are respectively L(y) defined by i i-> iy where xy is the Jordan 
product, and 

P(y) = 2(h(y)) 2 -IV) 

as defined in F.-K. page 32. If a and 6 are in V we denote by a ® b the endomorphism a; i— > atr(&x) 
of V. The endomorphism a ® 6 + b ® a belongs to J 7 . Denote by and T?, the linear subspaces of T 
respectively generated by and {d'y ® y + P(y); y € T^} and {y ® y — P(y);y G V}. From f|3 . 8() J 7 ! and 
Ti are canonically isomorphic to Q\ and Q2 by g i-> f q where q(x) = (f q (x),x). We endow T with the 
Euclidean structure defined by Tr(a&). Here again we distinguish the trace tr of the Jordan algebra V 
from the trace Tr of the endomorphisms on the linear space V. Here is a list of various traces: 

Proposition 5.1 

1. Tr(a®b) =tr(a&), Tr [(a ® b)(c <g> d)\ = tr (ad)tr {be), 

Tr ((ai (g) 61) • • • (a fc (g) 6 fe )) = tr (ai& fc )tr (a 2 6i) • • • tr (a fe fe fc _i). 

2. Tr [L(a)L(6)(c ® d)] = tr [(a(6c))d] 

3. Tr (P(a)(6 <g> c)) = tr [(P(o)6)c] 

Proof 1) is standard since it only involves the Euclidean structure of V and not its Jordan algebra 
structure. 2) is a consequence of 1). Applying the definition of P(a), 3) is a consequence of 2). 

In the theorem below, we consider an endomorphism of JF such that ^>(y ® y) = P(y) for all y E V. 
It is an essential tool of the two papers Casalis and Letac (1996) and Letac and Massam (1998). The 
theorem shows that T\ and Ti are its two eigenspaces and uses this fact to give the dimensions of the 
spaces of quadratic forms Q\ and Q2 defined in Th. 3.1 above. 

Theorem 5.2 

1. There exists a symmetric endomorphism ^ of J- such that ^(y^y) 

y{W{y)) = d'y®y+{l-d')f(y) 

2. The spaces T\ and T% are orthogonal and T = T\ © Ti 

3. The spaces T\ and T% are the two eigenspaces of \I/ corresponding to the two eigenvalues 1 and —d' 
respectively. 

4. The dimensions of T\ and JF2 are given by 

n(n + l) , , r(r-l) 1 + <f(2r - 3) + d' 2 (r - l)(r - 2) 

— — dim J 7 ! = dim J 7 ^ = — x £ ^ — 

2 2 1 + d' 

Examples. For the Jordan algebra associated to the Lorentz cone where r = 2 we get dim T<z = 1. More 
specifically, if E is a Euclidean space with scalar product x.y consider the Jordan algebra l^ = Mx£ 
endowed with the Jordan product between x = (xq,x) and y = (yo,y) defined by 

xy = {x yo + x.y, x y + y x). 

Here the Lorentz cone is {(xq,x) € V ; xq > ||x||}, the trace is tr(xo,x) = 2xq and the Peirce constant 
is d = dimE — 1. In this case J" 2 is spanned by the symmetry S defined by (x ,x) >-> (x , —x). To see 



= F(y) for all y € V. It satisfies 

(5.13) 
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this observe that if e = (1, 0) then S = e <E) e — P(e) is in T% and use dim — 1. As a consequence if 
a b 



b* 



c 



represents a symmetric endomorphism of V (where a is real, c is a symmetric endomorphism 



of E and b is a linear form on E) then 



a b 
b* c 

S = 



is in J-\ if and only if it is orthogonal to 



1 

-id B 



that is if and only if a = Trc. The dimension of J-\ is \{n — l)(n + 2). 

For the Jordan algebra S r of symmetric real matrices where d = 1, we get dim Ti — 15 ( r — l)( r + -0 
and dim J 7 ! = jr(r + l)(r 2 + 5r + 6). For the Jordan algebra of Hermitian matrices where d — 2, we get 

cBm 7- 1 = (^) 2 ,d te ,^ = ^- 1 » sl 

and since d' = 1, ^ is an orthogonal symmetry with respect to J-2. For the Jordan algebra of Hermitian 
quaternionic matrices where d = 4, we get dimJ^ = 4r r ( r ~ 1 )( r ^ 2 ) _|_ ^(r-i) anc j dimJ 7 ! = ^-(4r 2 — 1). 
For the Albert algebra where d = 8 and r — 3 we get dim = 27, dim J 7 ! ~ 351 = 27 x 13. 

Proof. 1) The existence of ^> is proved in Casalis and Letac (1996) (Lemma 6.3) and (|5.13[) is proved 
in Letac and Massam (1998) (Proposition 3.1). For proving that * is symmetric, enough is to see that 
Tr [ty(x <£> x) (y <g> y)] is symmetric in x and y in V since {y <£> y ; y € y} spans J 7 . Equivalently we have to 
see that Tr [V(x)(y <g> y)] is symmetric. From Proposition 5.1 part 3, we have to show that tr [(P(x)y)y)] 
is symmetric. Applying the definition of P, we get 

tr [{P(x)y)y)} = tr [{2{x{xy) - x 2 y)y]. 

Let us now use Proposition II. 1.1, (iii) in F.-K. which says 

L(x 2 y) - L,(x 2 )L,(y) = 2h(xy)h(x) - 2h(x)h(y)lL(x) . 

Applying this equality to y we get (x 2 y)y — x 2 y 2 = 2(xy) 2 — 2x(y(xy)) that we rewrite as 2(xy) 2 +x 2 y 2 = 
2x(y(xy)) + (x 2 y)y. Since the left hand side is symmetric in (a:, y) this proves 2x(y(xy)) + (x 2 y)y = 
2y(x(xy)) + (y 2 x)x which implies in turn that (2(x(xy) — x 2 y)y is symmetric in (x,y) and shows that VP 
is symmetric. 

2) and 3) Since {y®y; y G V} spans T and since 

y ® y = YTdf^'y ® v + p(y)) + TTd/^ ® y _ p(y)) ' 

clearly T = T\ + Ti- From the formula (15 . 13[) and the definition of ^ we get easily that T\ and J-2 
are made of eigenvectors of \& respectively for the eigenvalues 1 and —d'. In particular Jifl J2 = {0}. 
Therefore T = T\ © T2 and thus the endomorphism ^ has no other eigenvalues. From the fact that $ is 
symmetric, J 7 ! and J2 are orthogonal. 

4) It is the difficult point. We have dim J 7 ! + dim Ti — " ( '" 2 +1 ^ where n is the dimension of V. An 
other linear equation for (dim J 7 !, dim J^) is trace(^E') = dim T\ — d'dim Ti leading to 

dim T 2 = (^±12 tracc(v ^ . (5 . 14) 

We embark for a calculation of trace(\l/) by selecting an orthonormal basis / = (/^)™1™ +1 ^ 2 of J- and by 
computing Tr [^(ft)fe\ in order to get 

n(n+l)/2 

trace(*)= £ Tr [9 (ft) ft]. 
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The basis / is chosen as follows. We start from a Jordan frame (ci, . . . ,c r ) of V (see F.-K. page 44). 
Recall that c 2 = c s and c s Ct — for s ^ t. We denote by V(c, A) the eigenspace of V of L(c) for the 
eigenvalue A. For l<s<t<rwe denote 

v st = v{c si ~) n f fa, i), v ss = v(c„ l). 

Recall that V = ©i< s <+< r Vst, that the dimension of V s t is for s < t and 1 for s = t and that these 
spaces are orthogonal (F.-K. Th. IV 2.1, (i)). Let (c* t , . . . , cf t ) be an orthonormal basis of the space V s ,t 
for s < i. The space V ss is spanned by c s . For simplicity denote also by e = (ei, . . . , e n ) the orthonormal 
basis of V defined by the c' s s and the c* t 's. Finally the basis / of T consists of the elements of the form 
fi — e« <8> e, for % = 1, . . . , n, or = (e, §5 ej + ej ® ej)/v2 for 1 < £ < j < n. Since e is an orthonormal 
basis of the Euclidean space V it is standard to see that / is an orthonormal basis of the space J- of 
symmetric endomorphisms of V. 

We now compute Tr \^{fi)fi] = Ct for all possible choices of fi in the basis /. 

1. Case A: ft = e, ® e,. From Proposition 5.1, part 5 we have for all a; £ V : 

Tr(P(x)(x<g> x)) = trx 4 (5.15) 

Case Al: e, = c s . Thus inserting x = c s in (|5.15l) we get Ci = trc s = 1. 

Case A2: = cj t . We use the fact that x 2 = ^j-(c s + c±) when x £ V, t (see F.-K. Proposition IV. 
1.4 (i)) and apply (|5.15l) to x = c k st . We get 

Q = tr(4) 4 ^itr[( Ci+Cs ) 2 ] = i 

2. Case B: ft = (ei <£> ej + ej ® ei)/y/2. We use the following calculation: 

*(x ® y + y ® x) = P(x + y) - P(x) - P(y) = 2[L(x)L(y) + L(y)L(x) - L(xy)] 
(F.-K. page 32) and, using Proposition 5.1 part 2: 

Tr [(L(x)L(y) + L(y)L(x) - JL(xy))(x ® y + y ® x)] = tr [(yx 2 )y + (xy 2 )x]. 

Thus 

= tr [(e 3 e 2 ) ej + (e ie 2 )e 4 ]. (5.16) 
Case Bl: = c s , ej = c< with s < t. From (15.161) : 

C« = tr [(c t c 2 s )c t + (c s c 2 )c s ] = 0. 

Case B2: ej = c s , ej = with 1 < u < v < r with s £ {it, v}. By the definition of V uv we have 
c L c s = \ c uv an ^ thus from (j5.16|) : 

C t = tr [(c^c 2 )cL + {c s {c k uv f)c s ] = tr [(<*,<:.)<*, + (c s l(c u + c v ))c s ] 
= tr [-(c„ + c^) + (c s i(c u + c„))c a ] = 1. 

Case B3: ej = c s , ej = c k v with 1 < u < v < r with s g" {u, v}. Here we have c\ v c a = from F.-K. 
page 68 last line. A calculation similar to B2 gives Ce = 0. 

Case B4: e, = c\ v , ej = c™ with 1 < u < v < r and 1 < k < m < d. Note that if x and y have 
norm 1 in V uv then 

(yx 2 )y + (xy 2 )x = (y-(cu + c.„))y + (x-(c u + c„))x = i(y 2 + x 2 ) = -(c u + c v ) 
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Applying this to x = c\ v and y — c™ we get Cf — 1 through (|5. 16[) . 

Case B5: e, = c^ t , ej = c™ with 1 < s < t < r, with 1 < u < v < r, with 1 < k, m < d and with 
{s,t} H {u, v} reduced to one point, say u — s. Note that if x and y have norm 1 in x € U s t and 
y € V sv then 

(yx 2 )y + {xy 2 )x = (y^{c s + c t ))y + (x^{c s + c v ))x = i(y 2 + x 2 ) = i(2c s + c t + c v ) 
Applying this to x — & st and y — c™ we get Cg = 1/2 through (|5.16[) . 

Case B6: = c k stl ej — c™ with 1 < s < t < r and with 1 < u < v < r with {s,t} fl = 0. 

Using F.-K. page 68 last line we see that (yx 2 )y + (xy 2 )x — when x € V st and y £ V uv . Therefore 
C t = 0. 

We are now in position to compute the trace of 'J. We adopt the obvious notation C{Ai) — J2i£A ^i- 
Thus 

trace(*) = C{A X ) + C{A 2 ) + C{B 2 ) + C{B 4 ) + C(B 5 ). 

Since Cg is constant on each of these five sets Ax, A 2 , B 2 , -B4, B5 we hrst count the number of their 
elements: 

N(Ax) = r, N(A 2 ) = r(r - l)d', N(B 2 ) = 2r(r - l)d', 
N(B 4 ) = r(r - l)d'(d' - ~), N(B 5 ) = 2r(r - l)(r - 2)d' 2 . 

We get hnally 

traced) = r + r(r - l)d' [2 + (r - l)d'} 
which leads to the result through f|5.14|) . 

Comments. Observe that Tr [^(J)/] = Tr/ 2 if and only if / is in T\ (write f = fi + f 2 with fi e Ti 
and Tr [*(/)/] - Trf 2 = -{d! + l)Tr/| to see this). Thus in the above orthonormal basis (/£)™ ( ™ +1)/2 
of T we have fi e T\ if and only if Ct = 1, which happens only in the cases Al, B2 and B4. This is a 
set of size N(A t ) + N(B 2 ) + N(B 4 ) < dim T\ . Similarly Tr [*(/)/] = -d'Tr/ 2 if and only if / is in T 2 , 
and this shows that no fi is in T 2 since Cg > for all i. 
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